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ABSTRACT 

The origin of the prompt emission in gamma-ray bursts (GRBs) is still an unsolved 
problem and several different mechanisms have been suggested. Here we fit Fermi 
GRB data with a photospheric emission model which includes dissipation of the jet 
kinetic energy below the photosphere. The resulting spectra are dominated by Gomp- 
tonization and contain no significant contribution from synchrotron radiation. In or¬ 
der to fit to the data we span a physically motivated part of the model’s parameter 
space and create DREAM (Dissipation with Radiative Emission as A table Model), 
a table model for xspec. We show that this model can describe different kinds of 
GRB spectra, including GRB 090618, representing a typical Band function spectrum, 
and GRB 100724B, illustrating a double peaked spectrum, previously fitted with a 
Band-l-blackbody model, suggesting they originate from a similar scenario. We sug¬ 
gest that the main difference between these two types of bursts is the optical depth 
at the dissipation site. 

Key words: gamma-ray burst: general - radiation mechanisms: thermal - gamma-ray 
burst: individual: GRB 090618 - individual: GRB 100724B 


1 INTRODUCTION 


For decades, the prompt y-ray spectra of gamma-ray bursts 
GRBs) have been fitted with the empirical Band function 
Band et al.l 1 1991^ 1. Although often producing good fits to 
GRB spectra, the model does not represent any physical sce¬ 
nario. In order to extract physical information from the ob¬ 
servations we need to introduce physically motivated models 
and fit them to data. 

There are several reasons why such models should 
include emission from the jet photosphere: (i) Syn¬ 
chrotron radiation fails to explain the observed GRB spec¬ 
tra due to the line of death _|Preeceet_^ Il998l ) and 
spectral width llAxelsson &: Borgonovd 20151 ): (n) Some 
GRBs have s pectr a which are close to pure Planck func¬ 
tions ( Rvdel I2OO4I: iGhirlanda, Pescalli fc Ghisellinil l2013l : 
iLarsson, Racusin fc Burgess! l2015h : (in) Many GRBs are 
well described by models c omprising a blackbody (BB) and 
an additional component l|Rvde et al.l 20101: iGuiriec et al.l 
l201ll : lAxelsson et al.ll20lj : lBurgess et al. 20141 ). 

At the same time, most bursts produce spectra not con¬ 
sistent with the simplest version of photospheric emission; 
in general the low-energy slope is far too soft to be ac¬ 
counted for by a BB. However, it has been realised that we 


* E-mail: bjornah@kth.se 


should not expect a pure Planck function from the photo¬ 
sphere in a relativistic outflow, since there are natural mech¬ 
anisms whi ch broaden the spectrum, e.g. subphotospheric 


anisms wUi cU broaden tue spectrum, e.g. subpnotospUeric 
dissipation jRees fc MeszarosI 200^ Pe’er^ Meszaros fc Reed 


120061 : [Gianniosll2006l : Ghhotrav fc Lazzatill2015l). or g eomet- 

ric broadening I Pe’ei 20081 : Lundman, Pe’er fc Rvdel[20H ). 
In this work we consider a model based on the former and 
fit it to GRB data from the Fermi Gamma-ray Space Tele¬ 
scope. Subphotospheric dissipation has already been sug¬ 
gested as the emission mechanism in some bursts, e.g. for 


GRB 090902 B liRvde et al.l 1201 ll: iPe’er et al.l 120121) and 
GRB 110920 (llvAni et al. 20151) . but a full physical model 
has not previously been fitted to data. The analysis pre¬ 
sented in this Letter provides a first proof-of-concept of fit¬ 
ting such a model to data. 

Throughout this paper we are assuming standard A- 
GDM cosmology, a Hubble constant of Hq = 69.6 and a 
cosmological constant, Ha = 0.714. All errors are 1 sigma 
unless otherwise stated. 


2 THE MODEL 

The physical scenario we consider is based on the fireball 
m odel and the radiat i on is t reated using the code described 
in IPe’er fc Waxmaiil (l2005l ). For a review of the fireball 
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Figure 1. Spectra obtained from the code for different values 
of the r of the dissipation site. The other parameters are F = 
250, Lo, 52 = 10, £pi = 0, eb = 10“®,£e = 0.9 and = 0.2. In 
addition for each spectrum we also show the BB spectrum of 
the seed photons, plotted with the same style and colour as the 
spectrum, but with thinner lines. 


model, see e.g. iMeszaro i (|200filV In our picture a progen¬ 
itor releases a luminosity I/o ,52 = erg s“^ (not to 

be confused with the observed luminosity), in a relativis¬ 
tic, collimated and magnetised jet of electrons, baryons, and 
photons. The jet is accelerated up to the saturation radius 
rg such that rg ~ ro? 7 , where the bulk flow Lorentz factor 
T — r) and where t] = Lo,s 2 /Mc^ is the dimensionless en¬ 
tropy and ro the nozzle radius. We assume dissipation to 
occur at a radius Vd, defined by the corresponding optical 
depth T. A fraction edLo,52 of the energy is dissipated by 
some, in principle unsp ecified, process, e.g. inte rnal shocks 
(|Rees_fc.^Ieszar^|2005|) ^ magn etic reconnection l|ThomDsorJ 
I 1994 I: iGiannios fc Spruiul2005l i or hadronic collision shocks 
( Beloborodov! ^id i. Here we assume ro = thus us¬ 

ing i nternal shocks for this mechanism llMeszaros fc Re^ 
Il999li . In practice, the only effect of this assumption is that 
it sets the initial photon temperature, ro(ro). The dissipated 
energy is divided between magnetic fields, which receive a 
fraction ebedTo, 52 , and the electrons, receiving eeedTo, 52 . A 
fraction Epi of the electrons take on a power law distribution 
and a fraction (1 — £pi) assume a Maxwellian distribution. 

Photons and particles interact via Compton and inverse 
Compton scattering, pair production/annihilation and syn¬ 
chrotron self-absorption, and the electrons emit synchrotron 
radiation. The code follows the spectral evolution of the elec¬ 
trons and photons over one dynamical time, tdyn = rd/c with 
a fine time resolution. Although the code does not simulate 
the hydrodynamical evolution of the jet, we can evaluate the 
time evolution in a GRB by performing time-resolved spec¬ 
troscopy and assuming that the jet properties are driven 
by changes of the central engine. From the dissipation ra¬ 
dius to the photospheric radius there should be adiabatic 
expansion, effectively cooling the jet, however this effect is 
neglecte d since it is small and comparable to other uncer¬ 
tainties (jPe’er fc Waxmar]|2004l . 

In order to be able to fit this model to data we cre¬ 
ate a grid of models by running the code for different in¬ 


put parameters. The grid is then turned into a table model 
for XSPEC. We will refer to the model as DREAM {Dissipa¬ 
tion with Radiative Emission as A table Model). This way 
we may perform relatively fast fits even though the sim¬ 
ulations are computationally expensive. For this study we 
have chosen to confine ourself to a four-dimensional param¬ 
eter space, using the parameters r, F, Lo, 52 , £d, keeping the 
other three parameters fixed at eb = 10“®, £e = 0.9 and 
£pi = 0. This choice reflects a scenario where the vast major¬ 
ity of the energy goes to the electrons, which take on a com¬ 
plete ly Maxwellian distribution as they are heated llLevinsonI 
I 2012 !) . and where we have weak magnetic fields, yielding 
negligible synchrotron radiation. This allows us to address 
the question of whether the observed spectra can be ex¬ 
plained without this process. The values for the parameters 
used are r = 1,5,10,20,35, F = 50,100,250,500, Lo ,52 = 
0.1,1,10,100, 300, £d = 0.1, 0.2, 0.3, 0.4, 0.5, yielding a table 
model consisting of 500 simulations and spanning a physi¬ 
cally motivated part of the parameter space. In xspec the 
model obtains two additional parameters; a redshift, 2 and 
an additional normalization which is proportional to the ob¬ 
served photon flux. In order to eliminate degeneracies we 
keep these two parameters constant for each burst. Addi¬ 
tionally, to make sure that the resulting fits are not strongly 
affected by the step size in the table we also created a much 
finer grid spanning a smaller part of the parameter space. 
A comparison with the original model showed no significant 
differences. 

The shape of the resulting sp ectra for different input 
param eters has been discussed by jPe’er. Meszaros fc Reed 
(I 2 OO 6 II . Here we summarise the main points relevant in the 
context of our spectral fits. Note that the effects described 
below are non-linear and that the effect of changing one 
parameter partially depends on the values of the other pa¬ 
rameters. The code produces an output spectrum in terms 
of photon emissivity against photon energy in the comov¬ 
ing jet frame. In Fig. [1] we have plotted EFe spectra in 
observer frame for varying values of r. We note how the 
shape of the spectrum changes with increasing r, becoming 
more thermalized, ultimately approaching a Wien spectrum 
for higher values of r, as expected. Furthermore, we have in 
Fig. [T] included the initial BBs, corresponding to each spec¬ 
trum’s thermal seed photons. The figure illustrates that a 
low-energy spectral slope softer than Rayleigh-Jeans can be 
obtained as the thermal seed-photons are up-scattered. This 
effect is stronger for high optical depths. 

The luminosity, L 0 . 52 , corresponds to the amount 
of energy we have in the spectrum as well as the 
comoving proton number de nsity, since Lo ,52 oc rip 
ijPe’er, Meszaros fc Rees! l200d j. and thus a higher Lo ,52 
correspond s to a higher normalisation. C onsidering ra oc 
To.52/rr3 Ip e’er. Mdszaros fc Rees! [2006h . along with the 
photon temperature going as T = T'o(r’d/r'g)~^^^ 

llRees fc Meszaroal2005l ~l. we note that a higher Z/ 0,52 also 
moves the initial BB towards lower energies. In contrast, 
considering the dependence on F in the expression above, 
an increase in F results in fewer but more energetic photons, 
hence shifting the spectrum to higher energies and lowering 
the normalization. Another important effect of high F is pair 
production. When pairs are created in sufficient numbers 
they increase the effective optical depth, resulting in a con¬ 
siderably stronger thermalization. Lastly, £d increases the 
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temperature of the electrons, with the norm alised tempera¬ 
ture 6 oc Ed JPe’er, Meszaros fc Reesl[2006ll . yielding an in¬ 
creasingly higher second peak in the spectrum due to Comp- 
tonisation. Hence, as Ed approaches zero, the spectrum tends 
to the initial photon Planck function. 

Due to assumptions made in the code, we expect a gen¬ 
eral, underlying uncertainty in the parameter values that 
produce a given spectrum, estimated to amount up to a fac¬ 
tor of ~2. The main contributions to this uncertainty is that 
the code is 1 dimensional (and thus spatial effects are not 
considered) and that it does not a ccount for the hydrody- 
nami cal evolution of the jet (see alsol^ ’er fc Waxman|[20o3 . 
l2005h . It should be noted that this uncertainty is not cap¬ 
tured by the statistics when performing fits. 


3 OBSERVATIONS AND FITTING THE 

MODEL 

We illustrate the application of the DREAM model by fitting 
it to two different, bright Fermi GRBs; GRB 090618 and 
GRB 100724B. While GRB 090618 is well fitted by a Band 
function, GRB 100724B was the first ex ample of a burst wit h 
a significant additional BB component dGuiriec et al.l[201lll . 
GRB 090618 is analysed using G amma-ray Burst Monitor 
(GBM) data dMeesan et al.l 120091 ') from the Nal and BGO 
detectors. For GRB 100724B, we used GBM data from the 
Nal and BGO detectors as well as LAT-LLE data. The LLE 
data are a type of LAT data designed for the use with bright 
transients to bridge the energy range of the LAT and GBM 
dPelassa et al.ll2010l '). For both bursts we selected Nal de¬ 
tectors seeing the GRB at an off-axis angle lower than 60 
degrees and the BGO detector as being the best aligned of 
the two BGO detectors. The spectra were fitted in the energy 
ranges 8-1000 keV (Nal), 200-40 000 keV (BGO) and 30- 
1000 MeV (LAT-LLE). For a more extens ive analysis of the 
high e nergy data for GRB 100724B, see lAckermann et aP 
d2013l l who find evidence for an exponential cut-off at high 
energies. 

Our spectral analysis is time resolved, using a signal- 
to-noise-ratio (SNR) = 40 in order to bin the data. As a 
consistency check we have also performed the binning with 
Bayesian Blocks, and the results were not changed consid¬ 
erably. The analysis is performed using HEASARG’s xspbc 
12.8.Ig, with pgstat statistics and with PHA (Pulse Height 
Analyser) data. 


Table 1. Best-fitting parameters for the first five time bins for 
GRB 090618 fitted with the DREAM model. Cd is pegged at its 
lowest value. For the complete time evolution, the full table is 
available online. 


Time (s) 

pgstat/d.o.f. 

T 

r 

To, 52 

£d 

-2.0-4.0 

306/240 

14.0ti;« 

233 I 2 


0.100+°o°2 

4.0-8.0 

347/240 


242^3 

6.06t°;« 

0.100+°°°^ 

8.0-11.7 

294/240 

^•^^-0.41 

233/:3 

C 1^+0.29 

0.100+°°°2 

11.7-15.2 

247/240 

4.^U_Q gQ 

232tl 

O.UU_Q 28 

0.100+°°°2 

15.2-19.5 

304/240 

-0.27 

2A2+\ 

0 54+0.16 

0.100+°oo2 

Table 2. 

Best-fitting parameters for the first five 

time bins in 

GRB 100724B fitted with the DREAM model. For the complete 

time evolution, the full table is available online. 


Time (s) 

pgstat/d.o.f. 

T 

V 

To, 52 

Ed 

-1.0-9.4 

526/383 

n.btli 

4261^ 

69.8tti 

0.100+°°°2 

9.4-12.1 

454/383 


472t| 

1601" 

U.1ZU_Q QQg 

12.1-15.0 

406/383 

12.3Dd 

4691^ 

163t^ 

0.100+o°°2 

15.0-17.1 

481/383 

9.26i;®® 

470 I 5 

I 8 O/:® 

0 153+0015 

U.1 DO_q QJ 4 

17.1-19.4 

434/383 


475l| 

182t^ 

n 1 oQ+O.Oii 

U.100 Q Q08 


fit with a Band function, with the best-fitting parameters 
a = -0.92t°;°?, /3 = -2.33l°;12, Epe^k = 34011^ keV and 
the test statistic pgstat/d.o.f. = 236/240 0. 

3.2 GRB 100724B 

GRB 100724B has a fluence of 2.174 x IQ-"^ ± 5 x 
10“^^ erg cm“^, Tgo = 114.69 s llvon Kienlin et alll2014l l. 
and an unknown redshift. In order to conserve four free pa¬ 
rameters, we set the redshift to 1 when performing the fits 
to this burst. The most prominent effect of varying 2 for this 
burst is a change in Lo ,52 and £d due to the change in flux. 

The fit results for the first five time bins are presented in 
Tabled whereas the full table is available online. In Fig. 
we show a fit with the DREAM model at 25.8 — 33.5 s, 
after the trigger, as well as the corresponding fit with a 
Band-|-BB function. The DREAM model fit yields the pa¬ 
rameters T = 4 . 9 lio, 3 o, r = 44319 , Lo ,52 = 4I.9I2;? £d = 
0.121q qj with pgstat/d.o.f. = 406/383 and the fit with 
Band-l-BB yields the parameters a = —l.OOlQgg, /3 = 
-2.36lHt -Speak = keV, kT = 31.711,^ keV with 

pgstat/d.o.f. = 401/381. 


3.1 GRB 090618 


GRB 090618 has a fluence of 2.684 x 10“ ± 4 x 

10“^ erg cm~ ^ jvon Kienlin et al.ll20l3l and a redshift 2 = 
0.54 dCenko et al. I 2009 II , which translates into an observed 
luminosity L = 2.8 x 10®^ erg s“^. It is a long GRB with 
90 per cent of the flux emitted during Tgo = 112.39 s. The 
fit results for the first five time bins are presented in Ta¬ 
ble [U whereas the full table is available online. In Fig. [2] we 
present the bin at 65.3 — 65.7 s as an example. The fit with 
the DREAM model yields r = 17 .F = 239^2) -So ,52 = 
SS.OljJ £d = 0.100’^°'^'^^ and pgstat/degrees of freedom 
(d.o.f.) = 251/240. This test statistic shows that the 
DREAM model provides an acceptable fit to the data. For 
the same bin, also in Fig. [2] we present the corresponding 


4 DISCUSSION AND CONCLUSION 

As shown in section [3l we find that the DREAM model 
provides acceptable fits to the data. In addition, the pa¬ 
rameters are determined with small statistical errors, in the 
range of 3% — 23% for all parameters. These errors are of the 
same magnitude as the errors in the Band- or Band-|-BB- 
parameters. However, these errors must be assessed keeping 


^ Note that these test statistics should not be used for model 
to model comparison as the models are not nested. Furthermore, 
such comparison is not meaningful since only the DREAM model 
represents a physical scenario. 
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GRB090618 GRB100724B 

DREAM DREAM 



10 100 1000 10 ^ 


Energy (keV) 
GRB090618 
Band 




GRB100724B 

Band+BB 



Figure 2. Fit results for GRB 090618 and GRB 100724B at 65.3 — 65.7 s and 25.8 — 33.5 s, respectively. The first column shows 
GRB 090618; fitted with the DREAM model above and with the Band function below. Data from Nal 4 and BGO 0 are shown in black 
and magenta, respectively. The right column shows fits of GRB 100724B with the DREAM model above and with Band+BB below. 
Data from Nal 0,1 and BGO 0 are shown in black, red and magenta respectively. LAT-LLE data are shown in blue. The lower panels in 
these plots show the ratio between the model and the data points. Data have been re-binned for visual clarity, to errors of 4 sigma. 


in mind the much more essential uncertainties due to as¬ 
sumptions in the model, as noted above. In the case of ea 
it reaches its lowest allowed value of our examined param¬ 
eter space in most of the bins for these bursts. Hence, the 
corresponding one sided statistical errors, shown in Table [T] 
should be seen only as an indication of the fit’s sensitivity to 
an increasing ea. In our current configuration we find that 
no pairs of parameters are highly degenerate or correlated 
around the best-fitting values. Some degeneracies do exist, 
however, most notably between r and To, 52 , £a and Lo ,52 as 
well as r and r. The pegged value of ea implies that the value 
of To .52 might in reality be higher in these bursts than what 
is found in our fits. However, a significantly lower value of ea 
is not expected since that would yield a spectrum close to a 
Planck function, see the discussion in section [2] We have ex¬ 
amined progressively lower values of ca for the best-fitting 
parameter values of GRB 090618, which indicate that we 
should have ea > 0.01, as the observed spectrum is not close 
to a Planck function. However, due to the non-linearity of 
the problem, it is difficult to give a robust lower limit ea. 

When considering the fits in more detail we start by 
examining GRB 090618; this burst exhibits a typical Band 
function, in the sense that its best-fitting parameter val¬ 


ues lie close to th e average values of all GBM GRBs 
dGruber et aklEoidh . With a ~ —0.9 the soft part of the 
spectrum is certainly not what would normally be associ¬ 
ated with a BB. This is of particular interest as it is of¬ 
ten assumed that synchrotron emission is needed to obtain 
such values of a, whereas our model in its current state has 
no synchrotron component. Within the DREAM model the 
soft low-energy slope is instead created as a result of Gomp- 
tonization of the seed BB, c.f. Fig [T] In addition, in the 
Band fits the value of a depends on the assumption that 
the spectrum has a single peak. If the true spectral shape is 
instead weakly double-peaked (as suggested by the fits with 
the DREAM model, see Fig.m this approximation will push 
a towards softer values. 

GRB 100724B, in contrast to GRB 090618, is well fitted 
by a Band+BB model and thus represents a different cate¬ 
gory of burst. However, when fitted with the DREAM model 
both bursts appear slightly double peaked, see Fig. [5] This 
seems to suggest that these two bursts may originate from 
the same physical processes and should be considered to be¬ 
long to the same category. The reason why adding a BB re¬ 
sults in a greater statistical improvement for GRB 100724B 
than for GRB 090628 is likely driven by the distance be- 
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tween the two peaks, E 2 IE 1 , being significantly larger for 
GRB 100724B (qualitatively seen in Fig. [2]), with Ei and 
E 2 corresponding to the low- and high-energy peaks, respec¬ 
tively. Furthermore, the temperature of the seed BB deter¬ 
mines how much of the hard low-energy spectrum is seen in 
the GBM energy range, and thus whether a BB component 
can be identified (note that observed break in energy, Ei, 
does not correspond exactly to the seed photon energy), see 
Fig. [T] Within the DREAM model we find that, when ea 
is constant, the value of E 2 IE 1 is primarily set by r, with 
a higher t yielding a lower E 2 /E 1 . In this context we note 
that the optical depth may in reality be higher than what 
is given by r, due to pair production in the jet, which be¬ 
comes increasingly important with high F. This has no effect 
on any of our conclusions in this paper. 

An important conclusion that can be drawn is that, as¬ 
suming this model to be correct, the Band Flpeak does not 
represent the true peak energy of the spectrum, but is an 
average value of the two peaks. This point is something one 
would need to con sider when interpreting the Flpeak — E-iso 
llAmati et al.ll 200 ^ relation. Another feature of the model is 
the pair annihilation peak, seen at 10® keV for GRB 100724B 
in Fig. [21 Note that the placement of this peak is uncer¬ 
tain due to the unknown redshift. The DREAM model also 
predicts a sharp cut-off at higher energies, above the fitted 
energy range. 

In the evolution of best-fitting parameters of these two 
bursts there is a general trend that F varies little through¬ 
out the bursts, as one would expect from general fireball 
dynamics. The same is true for Ed, although we note that 
this parameter reaches its lowest allowed value in most of 
the bins, see discussion above. Furthermore the behaviour 
of Lq tends to approximately follow the light curves, as ex¬ 
pected from the nearly constant Ed- The evolution of r tends 
approximately to that of To, 52 , resulting in a weakly varying 
dissipation radius, see section[2] We note that the time scale 
for the parameter evolution is well beyond the time scale 
of the microphysica l processes of the emis sion. As pointed 
out by IPe’er) (l2008h : IPe’er fc Rvdd (l 201 ll 'l. there is an in¬ 
termediate time scale of light propagation in the scenario 
of subphotospheric dissipation, when taking e.g. geometric 
effects into account, which might be long enough to induce 
a temporal smearing between emission from multiple, dis¬ 
crete dissipation sites. This would imply a connection be¬ 
tween smoo th pulses in t h e ligh tcurve and subphotospheric 
dissipation. I Beloborodov! ll 201 lj 'l argues that this smearing 
is suppressed enough to allow for separate treatment of the 
individual pulses. In that case, as we assume that the dy¬ 
namical conditions are roughly constant during one pulse, 
the observed parameter evolution reflects the changing dy¬ 
namical properties of the central engine and/or dissipation 
process. 

In summary, we have shown that the DREAM model 
provides fits of acceptable quality to the data, using 
GRB 090618 and GRB 100724B as examples. Notably these 
fits has been performed with negligible synchrotron radi¬ 
ation. We have shown that these two bursts’ spectra are 
similar, having two spectral breaks as their main feature. 
We have thus demonstrated that they can be considered to 
originate from the same process, whereas previous fits to 
these bursts have claimed different components. The main 
difference between the two bursts is the distance between 


the peaks, which we chiefly attribute to the parameter r 
in our current setting. In a future work we will extend the 
model to include additional free parameters, in particular to 
include the Synchrotron regime. 
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